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We present a theoretical analysis of the nonlinear bias 
and temperature dependence of current-voltage characteris- 
tics of a spin-valve device which is formed by connecting a 
quantum dot to two ferromagnetic electrodes whose magnetic 
moments orient at an angle 6 with respect to each other. 
The theory is based on nonequilibrium Green's function ap- 
proach and focused on current perpendicular to plane geom- 
etry. Coulomb interaction has been taken into account ex- 
plicitly at the Hartree level. We derive a formula in closed 
form for current flowing through the device in general terms 
of bias and temperature. In the wideband limit we report 
exact results for the TMR junction nonlinear I-V curve as a 
function of 8. We also report the conductance slope at zero 
bias as a function of temperature for which experimental re- 
sults reported an anomalous behavior. 

73.23.Ad,73.40.Gk,72.10.Bg 



I. INTRODUCTION 

Due to the exciting perspective of magnetoelectron- 
ics where coherent charge and polarised spin are both 
utilized for electronic device functional, quantum trans- 
port properties of systems with a magnetic character 
have attracted extensive attention. The well known giant 
magnetoresistance effect (XiMR) is such a spin-polarized 
electronic transport effectBH, using which a GMR device 
changes its resistance depending on the orientation of 
magnetic moments of the magnetic multi-layers forming 
the device. A GMR system is fabricated by sandwiching 
non-magnetic metal layer between two magnetic layersH, 
and further extensions or variations of this structure can 
produce many other device functions, including the spift- 
valve transistoiQ, spin-selective electron interferometerO, 
and nonvolatile RAMEI. While GMR is perhaps one of 
the most fruitful areas of fundamental research which 
has already produced substantial application with huge 
commercial value, an entirely different approach, also ex- 
ploiting spin-polarized transport, is the tunneling mag- 
netoresistance (TMR) where an insulating (I) material 
layer is sandwiched in between two ferromagnetic (FM) 
layersQ, forming a FM/I/FM structure. One important 
attraction of a TMR device is that it carries lower current 
than the metallic GMR system, hence TMR can be, in 
principle, more advantageous for portable devices. Due 
to the reported room temperature operation of TMRE, 
the fundamental principle and various transport charac- 



teristics of TMR devices have been the subject of increas- 
ingly more research, and will be the topic of the present 
work. 

Tunneling of spin-polarized electrons through the 
TMR structure gives a large tunnelinSjHiagnetoresistance 
(TMR) through the spin- valve effectQiJ. The spin valve 
effect is characterized by the fractional change of resis- 
tance K = AR/R, where AR is the difference in resis- 
tance when magnetization of the two ferromagnetic met- 
als are in parallel or anti-parallel. In the original ex- 
periment of JuUiere on Co/Ge/ Fe, k was determined 
to be about 14% at 4.2K temperature at zero external 
bias voltage, and it dropped to 1% when bias was in- 
creased to-p. few meV. More recently, large TMR was 
reportedcHij to persist at room temperature where k is 
about 11.8% for CoFe/Al20^/Co at small bias. Again, 
in this dg"p<;e k decreases drastically when the voltage 
increasesQ't£l. These experimental results clearly point 
to the need of a thorough theoretical understanding of 
the bias as well as temperature dependence of TMR. 
On the theoretical side, TMR has been invpjtigated by 
many authors. Zhang et.al. have considered^ TMR us- 
ing a transfer matrix and calculated the tunnel conduc- 
tance and magnetoresistance of a FM/I/FM system with 
an oscillatory but decaying bias voltage dependence p£ 
TMR resulting from quantum resonance. Bratkovskyli3 
investigated a similar system and compared cases with 
or without impurity scattering. Barnas and Fert in- 
vestigatedptnagnetoresistance oscillations due to charg- 
ing effectalj; Zhang et.alX-^ investigated the zero-bias 
anomaly of the experimental dataEl and provided an ex- 
planation that the anomaly is due to localized spin exci- 
tations at|_t|h£.Jiamer-electrode interface; and a number 
of studiesM'E£ni3TlZl were devoted to improve the origi- 
nal TMR model of JulliereQ. Very recently, the nonlinear 
bias dependence of an interesting TMR device in the fopHi 
of FM/ 1 /FM/I/FM was investigated theoretically^. 
Due to the relative magnetic moment orientation of the 
three FM layers as well as the quantum well formed by 
the two insulating layers --sjarious very interesting TMR 
behavior can be expecteclE£l. 

The purpose of this work is to further investigate TMR 
effect as a function of external bias as well as temperature 
in a FM/I/FM device structure, where we pay special 
attention to the internal Coulomb potential due to inter- 
actions. For a quantum device subjected to an external 
bias voltage, the internal potential build-up can be an 
important factor determining transport characteristics at 



the nonlinear regime&HI, and this has not been included in 
previous analysis Ref. [ll|Jl^Jlq , |l4| of the bias dependence 
of TMR. There, external bias potentials are included non- 
self-consistently through Fermi functions of the electron 
reservoirs. In this paper we present a theory for TMR 
junctions based on the nonequilibrium Green's function, 
where the internal potential is solved self-consistently at 
the Hartree level. In particular we derive a general ex- 
pression for the electric current flowing through the TMR 
device in terms of external bias, temperature, and the 
angle 6 between the magnetic moments of the two FM 
electrodes. Due to the inclusion of the internal Coulomb 
potential into the theory, the predicted results are gauge 
invariant against the shift of the external bias — which 
is a fundamental requirement of any transport theory. In 
the wideband limit and applying the charging model for 
the internal potential, we further derive an analytic ex- 
pression of angular dependence of the I-V curve at zero 
temperature. The temperature dependence of TMR is 
then studied via a numerical integration of our analyti- 
cal formula. Finally, in the weakly nonlinear regime our 
theory allows predictions of the slope of the TMR near 
zero bias, where experiments^ showed a dip of magneto- 
conductance at low temperatures. Our results showed 
that the slope smoothly crosses over from a finite value 
at low temperature to zero at high temperature. 

The rest of this paper is organized as follows. In section 
II, we present the general theory by deriving a current ex- 
pression for TMR junctions. In section III, we apply this 
current expression for various situations, and a summary 
is given in section IV. 



In Eq. (|]), Hdot describes the quantum dot 

Hdot = ^endl^^dna ■ (4) 

Ht models the coupling between electrodes and the 
quantum dot region (the scattering region) with hopping 
matrix Tkcm ■ For simplicity we shall assume the hopping 
matrix to be independent of spin index, thus 
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In these expressions eka = e" + qVa with a = L,R; cj,^^ 
(with a =t, i or ±1 and a = —a) is the creation operator 
of electrons with spin index a inside the a-electrode. Sim- 
ilarly djjg. is the creation operator of electrons with spin a 
at energy level n for the quantum dot region. In writing 
down Eqs.(0) and (0), we have made a simplification that 
the value of molecular field M is the same for the two elec- 
trodes, thus the spin- valve effect is obtainedEI by varying 
the angle 9. Essentially, M mimics the difference of de: 
sity of states (DOS) between spin up and down electronsl 
in the electrodes. Finally, due to electron-electron inter- 
actions an internal potential build-up U{r) is induced 
inside the quantum dot region. Hence the actual Hamil- 
tonian of the quantum dot is Hdot + qU. 

To proceed, r-we first apply the following Bogliubov 
transformationE3 to diagonalize the Hamiltonian of the 
right electrode. 
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II. GENERAL THEORETICAL FORMULATION 

The TMR device we examine is schematically shown 
in the inset of Fig. (la). It consists of a quantum dot 
connected by two ferromagnetic electrodes to the outside 
world. The magnetic moment M of the left electrode is 
pointing to the z-direction, the electric current is flowing 
in the y-direction, while the moment of the right elec- 
trode is at an angle 9 to the z-axis in the x — z plane. 
In second quantized form this device is described by the 
following Hamiltonian, 



H = Hr 



Hr. 



H. 



dot 



Hi 



(1) 



where Hl and Hb, describe the left and right electrodes 
where a DC bias potential is applied. 
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4r. - cos(0/2)C+,, -as\n{9l2)Ctn, 
from which we obtain the effective Hamiltonian 
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The electric current flowing the TMR device modeled 
by Eqs. (|l|, R | 4 , |^) can be written, following the deriva- 
tions of RcfTpi and extending it to include spin depen- 
dent scattering, in terms of the Green's functions of the 
quantum dot (Ti = 1), 



Ic = iq fidE/TT)Tr {-2/m(S^) 
X [(G^-G°)/„ + G<]} 
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where /„ = f{E~ qVa) and the trace is over both 
the state index and spin index. Here G'^{E,U) is the 



2x2 matrix in spin space for the retarded Green's func- 
tion with C/(r) the electro-static potential build-up inside 
quantum dot region. In Hartree approximation G^{E, U) 
is given by 
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where the self-energy £»■ = S£(£: - qVL) + ^r[E - qVn) 
is also a 2 X 2 matrix in spin space, which describes cou- 
pling of the quantum dot region to the two magnetic 
electrodes. The detailed analysis of self-energy is given 
in the Appendix and the result is written as 
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with the rotational matrix Ra for electrode a defined as 

(13) 
= e and £!:„ 



r> _ f COS 6*^/2 sin 6*0/2 
' -sin6'„/2 cos6l„/2 



Here angle 6a is defined as 6'l = and 
is given by 
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The lessee. Green's function G^ in Eq. (|To|) is 
calculatedEil through the Keldysh equation Cj = 
G'~S<G'^ with the self energy S< = 1:<{E - qVl) + 
SJ^(£' — qVji) given by the following expression at equi- 
librium, 
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where Taa- = ~2/m(S^^) is the linewidth function. 

Using Keldysh equation and the explicit self-energy ex- 
pressions (Ol^Sl), after tedious but straightforward ma- 
nipulations one can confirm that Eq.(n(l) becomes 
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(16) 

Note that the same form of current expression is true for 
systems without spin dependent scatteringp3. But with 
such scattering, the trace is over states as well as spin 
variables. Eq. (|lq ) is the basis for our further analysis. 

Before presenting the results for the Green's function 
G"", it is worth to emphasize that it is now explicitly de- 
pendent on the internal potential landscape through the 
Hartree potential C/[(rX-||Itiis is an extension to the previ- 
ous NEGF analysisEyEScJ and it allows us to investigate 
TMR within gauge invariance requirement. At Hartree 
level U{y) is determined by the self-consistent Poisson 
equation 



V^C/ = A-Kiq I [dE /2tt) Y^{G<{E, U))c 



(17) 



Eq. (ITtI ) is a nonlinear equation because G'''° depends 
on U{r) by Eq. ([ll|) . Hence in this theory there is a need 
to self-consistently solve the coupled equations of G'''" 
and £/(r). This presents some analytical difficulties (see 
next section). 

Eqs. ( p^Jll| , p^ ) completely determine the nonlinear I- 
V characteristics of multi-probe TMR junctions. It is 
not difficult to directly prove that the current expression 
Eq.([l6|) is gauge invariant: shifting the potential at all 
the electrodes by a constant AV, Va ^ Va + AV, la 
from Eq. (|6|) remains the same. 

From now on we focus on the two-probe TMR device 
of the inset of Fig. (la) for which II = —Ir = I- After 
expanding the spin part of the trace in Eq. (|l6|), we 
obtain the final expression of electric current for the TMR 
device, 
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In this result, the 2x2 Green's function matrix G'"|-4S 
determined by the usual method of equation of motionEll. 
For our Hamiltonian the equation of motion can be solved 
exactly through straightforward algebra, and they are 
found to be: 
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A^A^-B^ 
where A^ and B are defined as 

A^ = E- Hdot ~qU-K 
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Before finishing this section, we note that Eqs. ( |lq , 
p^ , |2C| , ^ , p^ ) form the basis for numerical calculations 
of TMR I-V curves beyond what has been done before. 



since it is a gauge invariant formula. In a typical numeri- 
cal analysis, one computes Green's functions G,'' and the 
coupling matrix T using tight-binding modelsca; and the 
Poisson equation can be solved using very powerful nu- 
merical techniquesc3. Numerical analysis of these expres- 
sions is beyond the scope of this paper and in the next 
section we attempt to obtain the I-V characteristics via 
analytic means. 



III. APPLICATIONS 

As discussed in the last section, the spin dependent 
nonlinear I-V characteristics of the TMR device is ob- 
tained by self-consistently solving Eqs. ( p^p^ pO| , |2l| , 
|l7| ). Due to complexity of this problem, some reasonable 
approximations need to be made. In particular we will 
present applications of our theory from two approxima- 
tions. First, instead of solving the Poisson equation ( [l7| ) 
for C/(r) analytically which is perhaps impossible to do, 
we will parameterize U in terms of related geometrical ca- 
pacitance. This is the discrete potential, approximation 
commonly used in theoretical analysistjEj. This analy- 
sis allows us to calculate general nonlinear I-V curves at 
various temperature and qualitatively compare with the 
experimental observations. Second, at weakly nonlinear 
regime, the bias voltage is finite but small hence one can 
expand the entire current formula in terms of bias order 
by order. In this case the Poisson equation can be di- 
rectly solved to yield a spatially varying U{r). At weakly 
nonlinear regime, we calculate the second order nonlinear 
conductance which is..the slope of the I-V curve at zero 
bias. Experimentallya one observes a finite slope at low 
temperature, while zero slope at high temperature. 



A. TMR I-V curves and ratios 

To obtain a general qualitative behavior of the I- 
V curves, we shall simplify the current formula (llq) 
using the dis crste potential approximation. In this 
approximationEaO one assumes that C/(r) = C/q is a con- 
stant that provides a shift to the band bottom of the 
quantum dot due to Coulomb interactions. To solve for 
Uo, we introduce geometrical capacitance coefficients Ci 
and C2 between the electrodes and the quantum dot. The 
total charge AQ in the quantum dot is given by the right 
hand side of Eq. (O), and we parameterize AQ using Ci 
and C2, 
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Hence by evaluating the energy integration and using Ci , 
C2 as input parameters, one obtains Uo- 

The explicit expressions for the Green's functions and 
hence the current (ttq) and charge (G4), can be derived in 



the wideband limited where the coupling matrix P is inde- 
pendent of energy. In this limit the energy integrations 
of (18, |2J) can be completed. For instance, assuming 
T La = P/jct = Per , at zero temperature we obtain the I-V 
curve of the TMR device to be, 
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where P = P| + P^, AP = P| - P|, E^ = Ep - Eq - 
qUo + qVa, and 
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The internal potential Uo is determined by the following 
equation obtained by evaluating the charge AQ at the 
wideband limit, 
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To get physical insight we plot the I-V characteristics at 
different temperature parameter (3 and magnetic moment 
orientation 9. We have chosen P| = 1.0, P| = 0.4, Ci = 
C2 — 0.5 and Ep — Eq = —2 (units set by ?i = e = 2m — 

!)■ . 

Fig. (la) shows I-V curves at zero temperature /? = 00 
for different orientations 6' = (solid line), 0.47r (dot- 
ted line), and O.Stt (dashed line). Clearly current in- 
creases as the junction bias increases. When 9 = 0, i.e. 
when magnetic moment of the left and right electrodes 
are parallel, the current is largest at all bias voltages. 
When 9 = TT for which the moments are anti-parallel, 
the current is the smallest. This, is a well known result 
for TMR junctions at zero biasEl, but it holds at finite 
bias as well. This behavior can be seen more clearly in 
Fig. (lb) where we plot current versus angle 9 at differ- 
ent bias Ay = {Vi - V2) = 2, 5, 8, 11, 14. For all cases 
/ is minimum at 9 = n. When temperature is nonzero 
the energy integration of Eq. (llSf) cannot be completed 
analytically thus we integrate it numerically assuming 
phonons can be neglected. Fig. (2) plots current versus 
9 at different temperatures /3 = 1,0.5,0.2. The current 
is sensitively dependent on temperature and decreases 
substantially when temperature is increased (smaller /? 
corresponds to higher temperatureL|-|This is consistent 
with previous experimental findingsQH. 

To calculate tunneling magnetoresistance ratio, we use 
the definition of Ref. nS^ which defines a nonlinear re- 
sistance R = l/G = V/I and TMR ratio = (i?(0) - 



R{tt))/R{0). The TMR ratio versus voltage is presented 
in Fig. (3) where two curves corresponding to /3 = cxd 
and P — I are plotted respectively. At zero temper- 
ature, a large TMR ratio (about 30%) drops quickly 
as the voltage is-increased, a trend agrees with exper- 
imental findingsLrQ. At small voltages, the TMR ratio 
is much lower for higher temperature, also in agreement 
with experimental resulta3. We note that the theoretical 
values of the TMR ratio, about 34% at low lempera- 
ture, are much larger than the measured valuesE2l (about 
27% at 77K), due to our choice of ideal model and sys- 
tem parameters. But the qualitative trend of these re- 
sults are consistent. The TMR minimum in Fig. (3) 
is due to quantum resonance, as was also seen in previ- 
ous studiesEll. Finally, we note that experimental data 
of TMR as a function of bias is a slight concave-shaped 
curve peatjaero bias, but that of our model, as well as of 
othersE5lliHl3, are slightly convex. Thpi the decay rate of 
TMR is faster in the measuredpdatall3, which indicates 
that other physical mechanismscJ absent in the present 
model are perhaps needed to quantitatively understand 
the experimental data. 



B. Conductance dip at zero bias 

Experimental results on TMR junctions of Ref. @ 
showed that there is a clear conductance dip at zero bias 
at low temperature, a slight dip at intermediate temper- 
ature, and perhaps no dip at high temperature. This 
anomaly may be associated with a number of device de- 
tails such as the presence of impurities, localization, and 
scattering. In this section we will examine the zero-bias 
conductance dip as a function of temperature. 

The slope of conductance versus bias is the sec- 
ond order nonlinear conductance coefficient Gap-f = 
(Pla/dVpdV^ evaluated at zero bias voltages Vg — > 0, 
Ky — > 0. We apply a weakly nonlinear analysis of GafSj 
in which all quantities of interest are expanded order by 
order in bias voltage. This way one can directly solve the 
Poisson equation (|l7] ) without using the discrete poten- 
tial approximation: we seek the solution of U{r) in the 
following form. 
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where Ueq is the equilibrium potejitial and UQ,(r), Uap{v) 
are the characteristic potentialsESEa. It can be shown 
that jtiie-characteristic potential satisfy the following sum 



rulesHES ^^ Uq, = 1 and ^ o u^p = 0. Expanding G^ of 



Eq. ( |17D in power series of Va , from the Poisson equation 
(O) one can derive equations for all the characteristic 
potentials. In particular the expansions ace, facilitated 
by Dyson equation to the appropriate ordeiO: 
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with Gq the equilibrium retarded Green's function, i.e. 
when U = Ueq. At the lowest order we obtainEj 
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The first term on the right hand side of Eq. (p8|), which 
depends on internal potential Uq,, describes the induced 
charge density in the TMR junction. The second term of 
Eq. (Eq) is the injectivity which corresponds to the charge 
density due to external injection. Finally 
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and 
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are the local density of states. Once the characteristic 
potential is obtained from Eq. (P8|) , the the second order 
nonlinear spin dependent conductance Gap-f is given by 
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+ GSG;5/mS^GS)(/mS'-<5„^ - Imi::^){2up - 5^^)] (31) 

where up in Eq.(pl|) is a 2 x 2 diagonal matrix. 

The conductance slope at zero bias which has been 
measured in the experiments of Ref. y,|lO| is given by 
Gill, since Gm is the slope of dl/dV at zero bias. 
To simplify discussion we consider Gm near a resonant 
point where the Green's function is given by Gq ^ — 
l/iEp -Eo+ ir^/2) and G5,22 = 1/(^f - ^o + ^T^/^ 
In addition, apply the quasi-neutrality approximationEil 
we have ui = {dni / dE) / {dn / dE) . Thus expression (|3l| ) 
is easily calculated. In Fig. (4) we plot Gm as a func- 
tion of temperature at Ep — 2.5 while Eq = 3.0. As 
the temperature is increased {(3 decreases), Gm goes to 
zero smoothly. Our model thus indicates that conduc- 
tance slope Gm is nonzero at small temperature, and it 
diminishes to zero at. high temperatures, consistent with 
experimental resultsB. Our result also indicates that the 
crossover from a finite slope to zero slope as a function 
of temperature is smooth, as shown in Fig. (4). 



IV. SUMMARY 

In this paper, quantum transport properties of a meso- 
scopic conductor connected to two metallic ferromag- 
netic electrodes have been studied theoretically. We de- 
rived self-consistent equations which completely deter- 
mine nonlinear I-V characteristics of the TMR device. 
Our theory is gauge invariant due to the inclusion of 



long-range Coulomb potential. In the wideband limit 
the nonlinear I-V curves and the TMR ratio at different 
temperature and magnetic moment orientation 9 are cal- 
culated and our results are qualitatively consistent with 
the experimental measurements. We were also able to in- 
vestigate the conductance slopes at zero bias for various 
temperatures, where experimental measuijements showed 
an anomaly. Experimental measurementaj gave data for 
a few temperatures, and our result shows that it is a 
smooth crossover from a finite slope to zero slope as tem- 
perature is increased. Finally, in this paper we concen- 
trated on analytical predictions where several commonly 
used approximations were applied. Our results indicate 
that to completely explain all the anomalies in measured 
data of TMR, some further physical effects, which are 
neglected so far, may have to be included. 
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V. APPENDIX 



Using the Dyson equation again for the retarded Green's 
functions G^^ ^^ and G^^^ it' ^^ have the following rela- 
tions 

^fen,TT ^ 9kR]'^km,]]^mn,]T, ~^ 9kR]'^km,U^rnn,l] (36) 



The self-energy matrix Y7 has the following matrix 
elements from Eq.(|): (SOm'fe,TT = cos(0/2)T^*^„, 

sin(0/2)r* ,, (S'-W.ai = cos(0/2)r,V,. From 
Eqs.dl), (|3§), and (|3), we obtain 
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with the corresponding self-energy functions (E^)m„^|| 
and {^R)mn,u given by 

(Sfl)„„,|| = / ,Ti^RniTkRn[glmCOs'^{e/2) 



glj,^ Sin' (9/2)] 



(39) 



In this appendix, we will derive the self-energy SL^ 
For a = L, the self-energy can be written rather easilyEil 

[^L)mn,tTtT' — / , Tf.j^j-^TkLngkLa"a-a-' 

k 

where g^^^ is the free retarded Green's function for the 
lead a. However, for a — R, the situation becomes 
slightly complicated due to the spin-ffip process. In the 
following, wc will use-thc Dyson equation to calculate the 
self-energy functioned. To facilitate the calculation, we 
will define the following retarded Green's functions, 

G;„„,...(ii,t2) ^ -iO{h-t2){{d„,Ati),dt'(t2)}) (32) 



GL_,(ii,i2) ^ ^^e{h-t2){{CkRAtl),dt,{t2)}) (33) 

where m and n label the states in the quantum dot and k 
labels the states in the lead. To illustrate the calculation 
procedure, we consider GJ^„ ^ for simplicity. By using 
the Dyson equation 



— Urn 'T' UtqZj Lt 



(34) 



we obtain 






mm' ,\\^ra'k,\\^kn,\\ 



l<J-07mm',TT^m'fc,Ti^*:n,iT 



(35) 



(S5j)„„,ti = E nRmTkBm{glRX " gl.R^] ^0/2) COs(0/2) 
k 

(40) 
Similarly, we can obtain other self-energy functions 

(E^)„.„aT = E TkRmTkRnKRi ~ 91r^\ sHd/2) cos(0/2) 

k 

(41) 

(S5j)mn,ii = E '^kRmTkR7i[glR^ Sm^{9/2) 
k 

+ 5LiCOs2(0/2)] (42) 

Eqs.(|9|)-(||) are equivalent to Eq.(^. 
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fixed at w = 2. The other parameters are the same 
as those of Fig.(l). 

(3) The TMR ratio versus the voltage. Solid line f3 = 
cxD and dotted line (3 = 1. The other parameters 
are the same as those of Fig.(l). 

(4) The second order nonlinear conductance Gm, 
which is the slope of conductance at zero bias, 
versus temperature parameter /3. Other parame- 



ters: Fii = 1.0, Fix = 
Ef = 2.5, and Eq = 3. 
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FIGURE CAPTIONS 



Fig. (1) (a). The current / versus the voltage for different 
orientation angles: 9 = (solid line), 9 = OAtt 
(dotted line), 9 = O.Stt (dashed line). Inset: the 
schematic plot showing the TMR device considered 
in this work. (b). The current / versus orienta- 
tion 9 at different bias v = 2,5,8,11,14. Other 
parameters: F| = 1.0, F^ = 0.4, Ci = C2 = 0.5, 
Ep ~ Eq = —2., and h = e = 2m = 1. 

Fig. (2) The current / versus the orientation 9 at different 
temperatures: /3 = 1 (solid hue), /3 = 0.5 (dotted 
line) f3 = 0.2 (dashed Hne). Here the voltage is 
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